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FOREWORD 


This report summarizes the physical forces influencing Skylab 
experiments M551, M552 and M553 and was compiled under Contract 
NAS8-27015, This effort was performed by the Lockheed-Huntsville 
Research & Engineering Center for NASA -Mar shall Space Flight 
Center . 

The NASA Contracting Officer's Representative (COR) for 
this study is Mr. T. C. Bannister, S&E-SSL-T. The NASA princi- 
pal investigators for this study are Mr. R,M. Poorman (M551), 
S&E-ASTN-MM; Mr. J.R. Williams (M552), S&E-PT-M; and 
Mr. E.A. Hasemeyer (M553), S&E-PT-MWM. 

The author also acknowledges the other investigators whose 
results contributed to this report; C.M. Adams, University of 
Wisconsin; M.R. Brashears, Lockheed Missiles 8c Space Company; 
D.N. Braski, Union Carbide Corporation; J. L. Brown, Georgiy 
Institute of Technology; P.C. Johnson, A.D. Little Company; 

T.Z. Kattamis, University of Connecticut; D.J. Larson, Grumman 
Aerospace Corporation; K, Masubuchi, Massachusetts Institute of 
Technology; R.E. Monroe, Battelle Memorial Institute; J.M. Tobin, 
Westinghouse Electric Corporation; and W.A. Zisman, Naval 
Research Laboratory. 
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SECTION I. INTRODUCTION 


This summary report on physical forces is concerned specific- 
ally with the Metals Melting Experiment, the Exothermic Brazing 
Experiment, and the Sphere Forming Experiment. They were per- 
formed in the M512 facility during the Skylab I mission of May- June 
1973 by Astronaut Charles Conrad. These experiments are de- 
scribed briefly as follows: 

M551 Metals Melting Experiment : Three sample disks, each 
containing three metal specimens of varying thicknesses, are to be 
rotated automatically at a controlled speed of 2.5 rpm under an 
electron beam gun such that an electron beam weld seam is produced 
in the metal specimen at a radius of 6 cm. Disk materials include 
2219 aluminum, 321 stainless steel and tantalum. During the con- 
tinuous weld portion of each disk, both full and partial penetration of 
the disk will be achieved by having a constant power input but a vary- 
ing disk thickness. For each disk, the continuous weld will be 
followed by a dwell portion. In the dwell portion of the weld, the 
disk will remain stationary while the electron beam impinges on a 
thick segment of the disk, thus creating a large molten pool. The 
electron beam will then be shut off and the pool will be allowed to 
solidify . 

M552 Exothermic Brazing : A technique for joining stainless 
steel tubes will be tested in this experiment and the flow and solidi- 
fication behavior of weightless molten braze alloys will be studied. 
The joining technique will use a solid mixture that produces heat by 
exothermic chemical reaction to braze sleeves over 18.75-mm 
diameter tubes, using a copper -silver -lithium braze alloy. A pack- 
age containing four assemblies, each consisting of a tube with sleeve 
and preformed braze alloy surrounded by exothermic material, will 
be mounted in the M512 facility's vacuum chamber. The exothermic 
reactions in the four assemblies will be ignited in sequence and the 
whole package will be returned to earth for analysis. 

M553 Sphere Forming Experiment ; Twenty -eight 6. 3 5 -mm 
diameter spherical specimens will be cast using the electron beam 
gun as a heat source. The specimens will be initially supported on 
two wheels by a sting. After melting is completed, the spheres will 
then be separated from their stings and allowed to solidify while 
free-floating in the vacuum chamber. Specimens will consist of the 
following materials: pure nickel, Ni-1% Ag, Ni-30% Cu, and Ni-12% 
Sn. 



The M51Z facility consisted of a spherical processing chamber, 
approximately 40 cm in diameter in which several experiments were 
performed. The electron beam in the M512 unit was operated be- 
tween 50 to 80 mA at 20 kV. A vacuum environment of 10“'^ torr was 
created by venting the chamber directly to outer space. 

The nature, magnitude and direction of the physical forces 
which influence the experiments are important because they control 
the liquid dynamics of the molten metals. Fluid motions, in turn, 
affect the melting and solidification of the metals. Thus, the quality 
of the processed specimens depends on the physical forces and their 
attendant effects on liquid motions. 

For these particular Skylab experiments, the only significant 
difference between space and earth proces sing will be the lack of 
gravity. The maxirhum gravity level experienced during operation 
of the M551, M552 and M553 experiments aboard Skylab II was 
7 X 10" g^ (gE”9.8 m/sec^). Other environmental factors which 
may also differ from earth processing are the vacuum, radiation, 
electromagnetic and thermal conditions. 

Gravity has no direct effect on grain structure or other prop- 
erties of solidified material. These properties are determined by 
the crystallization kinetics which are controlled by short-range . 
intermolecular forces; i.e., the temperature and concentration at 
the fluid-solid interface. Gravity has not been shown to have any 
significant direct effect on these forces, but can affect solidification 
indirectly through its direct effect on fluid motion. The three major 
indirect effects of gravity on solidification are; 

• Sedimentation 

e Buoyancy-Induced Convection, and 

• Hydrostatic Pressure. 

In addition to these three very important indirect effects, the lack of 
gravity will also afford the study of a more direct effect — the oppor- 
tunity to obtain homogeneous nucleation. The long free-float times 
potentially available in space processing applications will allow 
melts to cool and nucleate without the deleterious effects of con- 
tainer walls. Wall effects are usually very strong in nucleation 
phenomena and prevent large degrees of supercooling from being 
attained in terrestrial processing. These effects are explained 
further below. 

Sedimentation : This effect may be significant whenever hetero- 
geneous mixtures exist in fluids, such as in monotectic, dispersed- 
particle, or fiber-reinforced composite casting. The denser of the 
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immiscible materials will tend to settle unless colloidal or electro- 
static attractions interfere. Also, in supercooled melts, segregation 
of freshly formed nuclei by gravity would affect the final grain struc- 
ture. Nonmetallic inclusions, gas bubbles and voids, which usually 
exist in melts, are also distributed nonuniformly by gravity. 

Convection: In terrestrial processing, gravity is the primary 
driving force for the convection of contained fluids when they are 
subjected to thermal or concentration gradients. Temperature gra- 
dients arise from external heating and cooling, whereas concentra- 
tion gradients are usually produced internally (Soret effect and solute 
rejection at interfaces which leads to "constitutional super-cooling"). 
These two gradients can produce density gradients large enough to 
induce buoyancy-driven flow. This fluid motion affects the tempera- 
ture and concentration profiles within the fluid. This subsequently 
alters the shape and rate of movement of the freezing interface, be- 
cause the kinetics of freezing depend on the local temperature and 
concentration. The degree of mixing caused by cpnvection may be 
large enough to change the rate of solidification from kinetic to heat 
transfer or diffusion controlled which could drastically alter the 
grain structure. Examples of this convective effect on the transitions 
from planar to cellular or dendritic growth, from columnar to equi- 
axed eutectic structures, and from unidirectional to colonied or 
banded eutectic structures are cited frequently in the literature. 
Another effect of gravity-driven convection on solidification processes 
is that the bulk fluid movement, if rapid enough, can break delicate 
dendrite arms and therby alter final grain structure. Furthermore, 
interlamellar spacing in eutectic growth and dendrite arm spacing 
are dependent on cooling rate which is a strong function of convection. 

Hydrostatic Pressure ; A body of fluid in a gravity field sustains 
a vertical pressure gradient as the bottom fluid must support the 
weight of the upper fluid. This pressure gradient distorts the shape 
of liquids on earth because the shape of a liquid surface is determined 
by the surface tension and the internal hydrostatic pressure (and ad- 
hesion if the liquid wets a solid surface). Distorted drops of liquid 
will result in nonsymmetrical solids upon freezing. 

Homogeneous Nucleation Effects ; Microgravity solidification in 
the M553 experiment may lead to containerles a freezing which offers 
the probability of homogeneous nucleation and attendant large degrees 
of undercooling . Most of the effects of large undercooling, and thus 
extremely rapid freezing rates, are beneficial. Increasing the 
amount of undercooling in an alloy leads to substantial reductions in 
both the segregation ratio of alloying elements and of dendrite-arm 
spacings. As a consequence, the rate of homogenization of an alloy 
increases with increasing undercooling. Other beneficial effects of 
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undercooling are reduction in grain size, refinement of endogenous 
inclusions such as silicates and sulfides, and homogeneous distri- 
bution of porosity. All these factors are beneficial to mechanical 
and possibly physical properties of an alloy, such as strength, 
hardness, electrical resistivity, etc. 

The following sections contain separate discussions for each 
experiment on the physical forces and, where applicable, their 
effects on molten metal flow and solidification, 
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SECTION n. M551 METALS MELTING EXPERIMENT 


This section contains analyses of the physical forces for the 
M551 experinnent. The first section contains an analysis of the 
driving forces, magnitude and pattern of natural convective fluid 
flows which occur in the Metals Melting Experiment. This section 
is followed by a discourse on the possibility of molten pool spattering 
due to the momentum of the impinging electron beam. Next the 
spiking in the partial penetration zone of the stainless steel samples 
is evaluated for gravity effects. This is followed by a discussion on 
beading phenomena. Finally expected dwell shape variations and 
measurements of surface forces for stainless steel are given. 

A. CONVECTION ANALYSIS 

As in the other M512 experiments, the physical forces and 
attendant fluid dynamics of the M551 experiment are the most impor- 
tant factors in determining quality of the final product (in this case, 
the solidified pool, seam and cut produced on various areas of thin 
metal disks). Flow patterns in the molten material are important in 
these experiments, because all of the materials have low entropies 
of fusion [1]. Thus, their solidification (microstructure) is con- 
trolled by the rate of heat transfer removal which changes with 

the Qow field [3 and 4]. The flow will be especially important in 
the dwell mode, since a relatively large pool of melt will be created. 
The degree of flow will also determine the amount of mixing attained. 
If no or little flow were present, all . heavier components would segre- 
gate to the bottom of weld zones on earth, but to a much lesser ex- 
tent in microgravity environments. Fluid flow can also affect the 
shape of the weld pool [7]. 

Application of dimensional analysis [5 and 6] to the govern- 
ing equations for eb welding, coupled with ground-based and KC-135 
experiments, should enable prediction of the extent of reduction or 
increase of motion in the weld pool and/or the change in flow pattern 
in electron beam welding in space. Possible physical forces which 
could induce fluid flow in the M551 experiment, and their causes, 
include:. 

• Effective Gravity Force : Resultant force on weld 
specimen due to earth's gravity and centrifugal and 
Coriolis forces of orbiting spacecraft, 

• Lor'entz Force : Electromagnetic forces induced by 
passage of the electron beam current through the 
specimen. 
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e Electrostriction ; Stresses induced when electrical 
permittivity changes with density. 

• Magnetostriction : Stresses induced when permeability 
changes with density. 

e Electrostatic Force: Caused by presence of excess elec- 
trical charge {due to beam current and/or thermionic 
emission), 

• Surface Tension ; Tangential stresses at vapor -liquid 
or liquid-liquid interfaces can be induced if surface 
tension depends on temperature and/or concentration. 

Surface tension will also cause pressure gradients 
across curved interfaces. 

« Density Differences Accompanying Phase Changes 

• Beam Force : Impinging electrons give up their 

momentum . 

e Thermal Expansion : Dilation and compression of fluids 
whose density changes appreciably with temperature can 
induce fluid flow. 

« Vibration : Uncontrolled movement due to engine opera- 
tion, astronaut motion, particle impacts , etc. 

• Centrifugal and Coriolis; Generated by disk rotation. 

© Vapor Pressure; Evaporating molecules impart 
momentum which leads to normal stresses at 
vapor -liquid interface, 

• Inertia Forces ; Tend to sustain induced motions. 

m Viscous Forces ; Tend to resist driving forces. 

The preceding forces, which could influence fluid flow and solidi- 
fication, appear explicitly in the conservation equations which apply to 
formation of a molten pool by electron beam heating. These equations 
are given in detail in Section IV. 

A formal method of determining the controlling physical forces 
affecting fluid flow and solidification in electron beam welding was 
introduced in Section IV. The controlling physical forces are deter- 
mined by nondimensioning the governing differential equations and 
performing an order -of-magnitude comparison on the various dimen- 
sionless groups which result. The key to successful analysis is in 
choosing the proper reference values; i.e., since no freestream velo- 
city exists, which forces do we equate to estimate a "characteristic" 
or "typical" velocity. Choosing the proper characteristic velocity is 
very important, since the reference time, temperature, etc., usually 
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depend on this velocity. The results arrived at in Section IV are also 
considered valid for the M551 experiment. 

The equation, in dimensionless form, governing electron beam 
melting is Equation IV. 1 shown in Section IV. Values of the pertinent 
dimensionless groups for each of the M551 materials are given in 
Table 1. Examining Table 1 in conjunction with Equation IV. 1, 


Table 1 


M551 DIMENSIONAL ANALYSIS 


Material 

^Oh 

"■bo 

^Oh 

P^o 1 

"oh 

''bo' 

^Oh 

2219 Aluminum 

5.8 X 10"^ 

3.09 

1713 

5287 

0.53 

321 Stainless Steel 

6.3 X 10'^ 

1.56 

1581 

2461 

0.25 

Tantalum 

3.4 X 10"^ 

7.19 

2942 

21150 

2.11 


*With Earth gravity (9.8m/sec^); ’‘'With 10"^ Earth gravity . 


an order -of-magnitude analysis indicates that surface tension driven 
convection will occur both in ground tests and for Skylab conditions 
because Nqh Furthermore, gravity driven convection will exist 

on ground tests, but will be negligible in the reduced gravity of Skylab. 
Thus different forces will control convection on earth versus Skylab. 
The preceding analysis also indicates that electromagnetic or Lorentz 
forces will be negligible with regard to causing fluid motion. See 
Section IV for further details. 

A review of the literature on electric arc welding, which is 
somewhat similar to electron beam welding, has generated the follow- 
ing facts which can be compared to the preceding conclusions. In a 
study of motion in weld pools in arc welding, Woods and Milner [7j 
conclude that Lorentz forces are the primary cause of motion. An 
additional secondary cause is the momentum imparted by the im- 
pinging arc. Surface tension forces were not considered. Kotecki 
et al. [8] showed that surface tension and momentum forces of the 
impinging arc controlled ripple formation in gas tungsten arc welds. 
Surface tension controls when the arc is shut off. Brimacombe and 
Weinberg [9] also conclude that surface tension is the driving force 
for fluid motion once the impinging jet is removed 
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One direct result of formulating the governing equations is that 
surface deflection (cutting action) is primarily caused by vapor pres- 
sure, but that the beam force is also appreciable. This was deter- 
mined from the vapor-liquid force equation of Section IV, which 
indicates that. 


Vapor Pressure 
Beam Pressure 


This agrees with earlier studies 

Examination of the 24 frames per second flight and ground 
movies for all three materials indicated agreement with the analyt- 
ical predictions presented earlier; 

© Significant molten metal motion was exhibited 
in both the dwell and continuous weld molds 
for both ground and flight tests. 

® Surface tension provided an equal magnitude 
of convection in the dwell pools in micro- 
gravity, but the filming speed was too slow 
to delineate flow patterns, 

B. WELD POOL SPLATTERING 

In both the M551 and M553 experiments, the heat released to 
the metal specimens by the impinging electron beam enables a molten 
pool to form. A question arises as to the stability of these molten 
liquids in low gravity. At least one of the KC-135 M553 specimens 
was seen to break up violently into many smaller liquid spheroids 
upon complete melting. The electron beam was still hitting the speci- 
men during break up. A similar instability might develop in the 
molten puddle formed during the dwell mode of the M551 experiment 
wherein the liquid might not adhere and separate from the solid disk* 
(violent splattering). 

There are at least three different mechanisms which might 
explain the instability observed in the M553 KC-135 specimens. Upon 
specimen may have experienced violent degassing which 
could have led to droplet breakup. Another mechanism might be 
electrohydrodynamic instability caused by interactions between the 
electromagnetic forces of the electron beam and the fluid flow field 
set up in the molten metal by both thermal gradients, Lorentz and 
other forces [lOj- Lastly, the momentum force associated with 
the impinging electron beam might have set up unstable surface 
oscillations on the molten metal. 


8 



The latter instability mechanism has been treated recently by 
Berghmans f 1 ij . He performed a theoretical study of fluid inter- 
face stability with special attention being given to the role of surface 
tension. It was motivated by its possible application to the splattering 
of molten metal as observed during electric arc welding. 

Berghmans' study concluded that the weld pool interface would 
be stable if the following condition was met, 

We^ < 1.04 + 3.3 Bo^ {II.l) 


where 


We = Weber number 
Bo = Bond number 


The analysis is only rigorous if inertia and viscous effects are small 
compared to surface tension effects. Inertia effects are not negligible 
in the M551 and M553 experiments, but the results of Berghmans' 
study should give a reasonable approximation. 

The criterion expressed by Equation (II.l) was applied to each 
of the materials in both the M551 (2219 aluminum, stainless steel, 
and tantalum) and the M 5 53 (pure nickel, Ni-Cu, Ni-Ag, and Ni-Sn) 
and for gravity levels between ground tests (gg = 980 cm/sec?) and 
those of Skylab (g = 10"^ Beam diameter was also varied be- 

tween 0.635 and 0.07 cm. For each of the above cases, 

2 -2 
We <10 

therefore, the momentum force of the electron beam will not be a 
primary cause of weld pool instability in either M551 or M553. 

Examination of the 24 frames per second flight and ground 
movies for all three materials indicated agreement with the analyt- 
ical predictions presented above. No unstable splattering was 
evident even in the microgravity environment of Skylab. 

C. GRAVITY EFFECTS ON SPIKING 

The electron beam welding phenomenon known as spiking occurs 
in the partial penetration welds of the stainless steel samples in the 
M551 experiment. This phenomenon is illustrated in Figure 1 
which indicates that spiking is caused by oscillation of the melt. 

Using a pinhole X-ray camera [l2] and radiographs of beam-on- 
plate welds [l3J to reveal beam-metal interactions in the cavity 
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FIGURE 1 


SPIKING MECHANISM IN ELECTRON BEAM 
WELDING 


foVsplLtogr” mechanic can ba propound 

With initiation of the beam, the beam vaporizes the 
material as it bores its way into the material. Full 
penetration is reached within 50 millisec with the estab- 
lishment of an equilibrium cavity. The beam, impinging 
on the cavity base, boils heated liquid up the walls of thf 
cavity, most likely in a whirling fashion. Complete metal 
mixing results due to this dynamic pumping action. Cavitv 
^osures are produced as the metal loses kinetic energy ^ 
Once a closure forms, it intersects the electron beam 
producing two effects. First, while the beam is heating 
the closure, regions below the closure are given time to 

b^ck i is heated and generally falls 

k into the cavity and the cycle recurs. However, the 
occasionally be vaporized in an explosion, 
the explosion, the cavity is clear of liquid and 

a snike”^ H again free to impinge on the cavity base causing 
back^nto thrc^^-f’ niore frequently the closure falls 

a tvDical rate iTri cyclic mechanism continues (at 

process! ^ equilibrium 

of is a failure or interruption 

fails to mittr^- mechanism. When a closure 

fri^P it ®^Pioded out, the electron beam is 

free to dwell on the cavity base for a longer time thus pro- 
ducing a penetration spike." fl2j ^ 

tra!L^ formed each time the electron beam pene- 

o he base of the weld and the severity of the spiking 
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depends on the frequency of the oscillation and the 
welding speed as well as the material being welded. 

"When the molten metal flows into the cavity, it 
can easily trap bubbles of contaminant gases at the base 
of the weld and, if the cooling rate is sufficiently rapid, 
the fluid may freeze before the gases can be convected 
to the surface, thus giving rise to porosity in the fusion 
zone. Similarly, if the cooling rate is sufficiently rapid, 
then the melt, as it falls into the cavity, may interrupt 
the beam long enough for the walls of the spike to freeze 
so that the fluid, upon falling into the spike, might not 
bond, thus forming a cold shut or crack. 

"The alternate penetration and closing of the cavity 
are also very efficient in mixing the meU and results 
in a very homogeneous fusion zone." [l3j 

A theoretical analysis of the melt-beam dynamics described 
above and in Figure 1 shows that spiking frequency is predicted 
by Reference [l3j. 



where 

h = maximum depth of penetration 
2a = weld width 
g = gravity level 

Py = vapor pressure of molten metal 
a = surface tension of molten metal 
U) = spiking frequency 


(II. 2) 


The preceding analysis considered welding a horizontal plate from 
above. The configuration utilized for M551 ground tests (beam im- 
pinges vertical plate from the side at a nine o'clock position with 
clockwise rotation of the plate past the immobile beam) also lends 
itself to cavity oscillations as the trailing melt tends to be above the 
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beam. It should also be pointed out that the above oscillation fre- 
quency represents the natural frequency of the cavity. A lower 
frequency, superimposed over the natural frequency, may occur 
due to an interaction of the static cavity closing forces (hydrostatic 
head and surface tension) and the dynamic forces. The dynamic 
forces arise from the motion of the wave or ripple which is formed 
in the melt as the cavity penetrates. 

Thus Equation (11.2) should not be expected to generate exact, 
rigorous results to match actual M551 experimental data. It can, 
however, be utilized to predict a low gravity variation in the stainless 
steel electron beam welds. As predicted in Equation (11.2), spiking 
frequency varies inversely with gravity level. Using the experimental 
conditions, physical properties and beam parameters for the M551 
stainless steel experiment and assuming 11Z0°C of superheat in the 
melt. Equation (II. 2) predicts oscillation frequencies of 22 Hz and 28 
Hz for gravity levels of 9.8 m/sec^ and 9.8 x lO”'* m/sec^, respec- 
tively. At a disk rotation rate of 1.58 x 10“2 m/sec at the beam 
impingement point, the spikes would be located 0,72 mm apart and 
0.56 mm apart in gravity levels of 9.8 m/sec^ and 9.8 10"'* m/sec^, 
respectively. This gravity effect is rapidly dissipated at higher 
superheats, however, as at 1270°C the spiking frequencies are 41 
Hz and 44 Hz at 9.8 m/sec2 and 9.8 10"'^ m/^ec2, respectively. This 
expected gravity-related variation in spiking will be very difficult to 
observe because it is very difficult to section radial welds precisely 
along their centerline and even more difficult to etch, polish and per- 
form metallographic analyses on curved sections. No evidence of this 
possible variation has been reported to date in any investigations being 
conducted by other contractors or NASA. 

D. BEADING CONSIDERATIONS 

Upon examination of the M551 disks, it has been noted that a 
beading occurs behind the electron beam as the beam moves through 
the molten material. At present, no detailed analysis has been con- 
ducted to explain this phenomenon; however, some consideration hs-s 
been given to understanding the mechanism. In order to present a 
hypothesis explaining the proposed mechanism some background 
information is presented. 

The movement of the electron beam cavity through the molten 
metal resembles the motion of an infinite cylinder (especially if full 
penetration occurs) through a fluid. It has long been known that such 
a body leaves in its wake a regular pattern of vortices which move 
alternately clockwise and counterclockwise and is known as a Karman 
vortex street. A schematic is presented in Figure 2, 
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FIGURE 2. KARMAN VORTEX STREET (DIAGRAMMATIC); 

STREAMLINES DRAWN IN A SYSTEM OF CO- 
ORDINATES MOVING WITH THE VORTEX 
STREET 


There exists a distinct frequency at which the vortices on one 
side of the wake are shed that deperida only on the Reynolds number 
of the fluid motion and is given by the Strouhal number, i.e., 


where f is the frequency, D is the diameter of the cavity and v is the 
velocity of the beam in the liquid metal. 

From experiment, the relationship of the Strouhal number to the 
Reynolds number has been determined and is shown in Figure 3. This 
shows that a critical Reynolds number exists below which no shedding 
has been observed. 

Using these concepts the theoretical number of beads per unit 
length has been determined for the M551 materials. It should be 
mentioned that no attempt has been made to account for the solid 
material (edges) in the cut or the extent of melt in the radial direction 
of the disks (which is a function of the thermal conductivity). These 
results are shown in Table 2, The surface tension is included as the 
magnitude indicates the tendency of the vortices to bead before solidi- 
fication occurs. 
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FIGURE 3. STROUHAL NUMBER VS 
REYNOLDS NUMBER 


TABLE 2 

BEAD POPULATION 


M551 

Materials 

R eynolds 
No. 

Strouhal 

No. 

Shedding 
F req. 
(beads/ 
sec ) 

Beads 
per cm 
(on each 
side ) 

Total 
Beads 
(both 
sides ) 

Surface 
T ension 
(dynes/ 
cm) 

321 Stainless 
Steel 

55 

0.125 

0.8 

0.5 

1.0 

1750 

2219 Aluminum 

86 

0.155 

1.0 

0.6 

1.2 

737 

Tantalum 

207 

0.19 

1.2 

0.75 

1.5 

2150 


Some qualitative observations can be made concerning these re- 
sults, First both aluminum and stainless possess Reynolds numbers 
near the critical values; thus it would not be surprising if some of the 
effects not accounted for could have a significant effect on the popula- 
tion of the beads if not controlling whether or not beading will occur. 
Tantalum, however, possesses the largest Reynolds number and value 
of surface tension; thus it is expected that these properties produce the 
greatest probability of beading and will give rise to the greatest popu- 
lation of beads. ^ 
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Average bead spacing distances were estimated for both flight 
test and ground based test specimens by inspection of photographs. 
These experimental values are compared to the predicted values 
(Table 2) in Table 3. 


TABLE 3 
BEAD SPACING 


M551 

Materials 

Bead Spacing, cm 

Theoretical 

Flight 

Test 

Groiind Based 
Test 

321 Stainless Steel 

1.0 

1.1 

1.1 

2219 Aluminum 

0.8 

— 

— 

Tantalum 

0.7 

3.6 

2.2 


The experimental bead spacing on the stainless steel specimens is 
seen to be in remarkably good agreement with the predicted value for 
both flight and ground based tests. No beading is exhibited on the 
aluminum specimen, however, and the spacing on the tantalum speci- 
men is considerably larger than the predicted value. Consideration 
was given to the possibility that the departure from vortex shedding 
theory might be due to "wall" effects resulting from the limitation 
on the lateral extent of melting. The presence of confining parallel 
walls will tend to increase the drag force on the cylinder due to the 
increased fluid velocity in the space between cylinder and wall. This 
may be considered equivalent to an increased "effective" Reynolds 
number for drag on the cylinder. Hypothesizing that the Strouhal 
number relates to this effective Reynolds number according to the 
relationship in Figure 3, the vortex shedding frequency should be 
increased by the presence of the confining walls. The results in 
Table 3 , however, indicate the opposite trend. Wall effects, there- 
fore, do not account for the observed departure from vortex shedding 
theory. Also, nearly identical bead spacing is noted for both flight 
and ground-based tests with the stainless steel specimens, while 
spacing on the tantalum specimens is larger for flight than for ground- 
based test. It is believed that the mechanism of vortex shedding is 
thus related to the beading pattern [l4J. The discrepancies for the 
aluminum and tantalum are presently under further investigation. 
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DWELL SHAPE 


E . 


Lack of hyrdrostatic pressure and dominance of surface tension 
will allow dwell pool shape to be different (more spherical) and mav 
cause shallower eb penetration. The different pool shape will occur 
because surface tension acts to pull the melt toward cooler surfaces 
which is away from the electron beam in this case. Assuming that the 
net gravitational field is always aligned with the beam in Skylab proces 
sing, the situations illustrated in Figure 4 will result regarding bulk 
melt behavior of the dwell pool. As seen in Figure 4a surface ten- 
sion dominance will pull the melt away from the beam, resulting in 
less resistance and deeper beam penetration than in terrestrial proces 
Sing when heating from above. When heating from the side, however 
terrestrial processing allows the melt to run out of the pool which will 
probably lead to even less resistance and more penetration than in 
space processing. Shallower dwell penetrations should be expected 
aboard Skylab because ground samples were processed as in Figure 


. As predicted, the gravity sag illustrated in Figure 4 was evi- 
dent in the ground samples shown in Figures 5a and 6a which also 
lead to more extensive dwell penetration than in the flight samples 
Wo sag was evident in the aluminum and steel flight samples shown 
in Figures 5b and 6b as the dwell shape appeared similar to Figure 


F. SURFACE TENSION MEASUREMENTS 


surface tension of stainless steel was experimentally 
determined (Ref. 15). The resulting values for the estimated 
iiquid/vapor surface energy are given below. 


T emperature 
(°K) 

■ 

Surface Energy 

(10'^ N/m) 

Liquid Density 
(gm/cc) 

1868 

1800 

6.97 

1943 

1779 

6.89 

1993 

1756 

6.81 

2058 

1710 

6.73 


adequate for engineering evaluation analysis. The 
h be estimated, and no previous measurements 

cfent o7vaSa? ^he important temperature coeffi- 

0 47^v?n-5 xiT surface energy for SS 304 aUoy is rather large, 

-0.47 (x 10 ^ N/m'^K) (erg/cm^oc). ^ ' 
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Microgravity 



a. During Space Flight 



c, During Ground Test if Heating from Above 


FIGURE 4. PROPOSED M551 DWELL POOL BEHAVIOR [l7] 
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a. Ground (Gravity downward) b. Flight (Microgravity) 

Fig. 5 Cross section of aluminum dwell region 



a. Ground (Gravity downward) 


b. Flight (Microgravity) 


Fig. 6 Cross section of stainless steel dwell region 
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G. 


CONCLUSIONS 


The primary conclusion which can be reached from the results 
now available is that no significant practical differences exist between 
terrestrial and microgravity electron beam melting. The differences 
found to date — oxide free front surface for the full penetration weld in 
the aluminum flight samples and the lack of gravity sag and shallower 
penetration in the stainless steel flight sample dwell — do not seem to 
have any practical benefits or detriments for conducting this operation 
in space. Although no results have been reported to date, the change 
in spiking frequency with gravity level and degree of melt superheat 
could lead to an increase in porosity and cold shuts in low gravity 
electron beam welding. 

It is noted that the only other known low gravity electron beam 
welding experiments ever conducted [l6j have come to similar con- 
clusions using similar materials, beam parameters, chamber and 
disks, movie coverage and low-g aircraft experimentation. Their 
conclusions [l 6] were that not much difference existed between Earth 
and low-g specimens; aluminum rendered continuous beading, where- 
as stainless steel resulted in droplet beading; weld shape and penetra- 
tion were the same on Earth and in low-g; and more porosity existed 
in one of.the low-g specimens. 

No qualitative evaluations could be made on the magnitude and 
pattern of flow in the dwells because filming speed was too low. This 
was true for both ground and flight tests. 

Surface tension and vaporization forces were dominant for this 
experiment both on the ground and in Skylab. Gravity forces were 
also important in determining dwell shapes for stainless steel and 
aluminum ground tests. 
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This section contains descriptions of the physical forces and 
their effect for the M552 experiment. The first portion of the chap- 
ter addresses the capillary flow analysis, while the second segment 
discvisses bubble dynamics associated with ^/1552, The next two sec- 
tions discuss variable gap flow effects and pertinent surface tension 
measurements, respectively. Finally discussions of flight results 
are presented 

The solid materials being brazed were tubular sections of 304L 
stainless steel and pure nickel. The brazing alloy was a eutectic 
composition mixture of silver and copper (72% silver, 28% copper) 
containing a small concentration (0.5%) of lithium. There were four 
specimen designs in which the most important variations were in the 
composition of the solid (stainless or nickel) and dimensions of the 
capillary gap between the inner tubular section and the concentric 
outer sleeve. In all specimens the silver-copper braze alloy was pre- 
placed in ring-grooves situated for that purpose, one at each end of 
the specimen. Upon melting, flow of the liquid silver-copper alloy 
into the gap was under the impulsion of capillary action, in which the 
force producing flow is primarily surface tension. Gravity, when 
present, has a modifying effect on the surface tension force field 
which may or may not be negligible, depending on the magnitude of 
the capillary gap. The four specimen designs are listed below; 

0,005" - Gap - 304L Material 

SLSl - Sky lab Flight Specimen 
MCS-1, 2 and 3 - Ground Characterization 
Specimens 

0.010" " Gap » Pure Nickel 

SLN2 - Skylab Flight Specimen 
MCN-1,2 and 3 - Ground Characterization 
Specimens 

0.020" - Gap - 304L Material 

SLS3 Skylab Flight Specimen 
MCS-4, 5 and 6 - Grovmd Characterization 
Specimens 

0.000" - 0.030" Taper Gap - Pure Nickel Material 

SLN4 - Skylab Specimen 

MCN-4, 5 and 6 - Ground Characterization 
Specimens 
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A. CAPILLARY FLOW 

The effect of gap width and low gravity on capillary flow in the 
M552 Exothermic Brazing Experiment has been studied. Time-to- 
spread, flow velocities and Reynolds number have been calculated 
for the configuration shown in Figure 7. The results are shown 
in Tables 4, 5 and 6. 

The results indicate that [l7j; 

• Spreading time will be up to 5 0% shorter in gravity 
fields anticipated for Skylab than in ground tests. 

e Turbulent or oscillatory laminar flow will occur in 
Skylab processing in certain gaps where only laminar 
flow could occur in ground tests. 

The flow equation for an annulus is given by, 


where 

= velocity 

= axial position in annulus 
= density 
= sleeve radius 
= viscosity 

= tube radius/ sleeve radius 

_ 1 - k^ . 1 - 

' 1 . 

g = gravity 

V, - 2a cos 9 
o ~ p g R (1 - k) 

a - surface tension 
9 = contact angle 

Equation (III. i) was developed assuming laminar, steady state, incom- 
pressible, isothermal flow with g opposing the flow direction. It 


u 

h 

P 

R 

P 

k 

0 
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le 



Slot for 
Braze Ring 


FIGURE 7. M552 CONFIGURATION 


Table 4 

M552 SPREADING TIMES 


Gravity 

Gap W idth 

Spread Time (Annulus) 

(k/Ke;) 

(in.) 

{jl sec) 

2 • lO’^ 

.005 

12.00 


.010 

6.00 


.020 

2.97 


.030* 

1.76 

1 

.005 

13.60 


.010 

lAl 


•020* 

4.42 


.030 

3.60 


Assuming no variation in gap width. 

Table 5 

M552 FLOW VELOCITIES 


Geometry 

Gravity 

(g/g£) 

Annulus 

2 • 10“^ 

1 

1 

Slot 

2 • 10"^ 


1 


Gap 

Un.) 


.005 

.010 

.020 


Velocity (u) 
(cm/ sec) 


Assuming no variation in gap width 



























Table 6 

M552 REYNOLDS NUMBERS 


Geometry 

Gravity 
(g/ Be) 

Gap 

{in.) 

Rei 

Flow 

Regime 

Annulus 

2-10*^ 

.005 

.010 

.020 

.030^ 

202 

794 

3200 

7210 

Laminar 

Oscillatory 

Turbulent 

Turbulent 

1 

.005 

.010 

.020 

.030* 

182 

650 

2170 

3490 

Laminar 

Laminar 

Oscillatory 

Turbulent 

Slot 

2.10"^ 

.005 

.010 

.020 

.030* 

99 

395 

1580 

3550 

Laminar 

Laminar 

Laminar 

Turbulent 

1 

.005 

.010 

.020^ 

.030* 

90 

315 

970 

1500 

Laminar 

Laminar 

Laminar 

Laminar 

. 


Assuming no variation in gap width. 




should be a valid approximation for the M552 configuration and pro- 
cessing conditions. 

For very narrow gaps, Equation (III.l) reduces to that for flow 
in a slot which is given by 


u 


^ g 
12m h 



- h) 


where 


D = slot width 

This equation was developed for assumptions similar to those of 
Equation (111.1) and will be just as rigorous except for geometrical 
effects. 

The Spreading time, t, for flow in an annulus can be given by 


t (h„£n(l -h/h^)+h) 

pgR (j) 

Again, thi^ equation is valid for isothermal, laminar, incompressible 
flow with g opposing flgw. Levich [l8j states it will only be an ap- 
proximation if t < — . 

The Reynolds number for flow in an annulus is given as Reference 


Re = up 

1 M 

or by Reference 20, 

4R^ up 


where R^^ = hydraulic radius. 

The transition from laminar to turbulent flow occurs at Rco 2300. 
The transition from laminar viscous flow to oscillatory flow is given 
by Re 2 = 700. By making the transformation, k = 1 - e , and allowing 
e — *>-0, it can be shown that Re^ — ► Re 2 as k — ► 1. Thus the 
transition points for Rcj^ are approximately 700 and 2300 also. The 
Reynolds nximber for flow in a slot is defined as 
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Re 


s 


p u D 


and the laminar -turbulent transition is at Reg w 2000. 
B. BUBBLE DYNAMICS 


Many ground and KC-135 aircraft tests for M552 have exhibited 
a bubble or void space within the braze structure. Analysis of the 
movement of such a bubble in low gravity and the nonisothermal con- 
ditions of M552 has been completed. These analyses are applicable to 
motion of gas bubbles only and will not predict the location of shrink- 
age voids. 

Neglecting wall effects, the limiting velocity of a single bubble 
moving in a homogeneous liquid \mder the combined influence of a 
gravity field and a temperature gradient is given by Reference 21, 


2 Tt a da dT , , 2 , . , 

“ = T7TB- P dT dT- - <'’g - Pi!> 8^ ('‘g'*-"P 

8 *■ 

where 

A = 2 + k,/k 

r g 

B = 3Atg + 2fi^ 
k = thermal conductivity 
/u = viscosity 
a = bubble radius 
a = surface tension 
T = temperature 
p = density 
g = gravity 

and subscripts f and g refer to the liquid and gas phases, respectively. 
The temperature gradient affects bubble motion through tangential 
stresses at the bubble-liquid interface. These stresses are caused by 
surface tension gradients across the bubble. The surface tension 
mechanism is termed the Marangoni effect, which can be caused by 
either temperature or concentration gradients. 
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The ratio of buoyant /Marangoni forces on the bubble is given 


by, 


(p^ - P^) ga (Pg + pp A 

V = — — V"li 

^ ^ g dT dz 


where y = buoyancy/surface tension forces. Values of y for various 
size bubbles and gravity levels are given below for M552 conditions 
using a 10°C/cm temperature gradient, 


Bubble 



Diameter 

Gravity 

r 

(cm) 

(g/gE> 


10 -" 

1 

10 ^ 

10 “^ 

-6 

10 

1 

10 "^ 

1 

10 ^ 

10-2 

10"^ 

10”^ 


Thus, Marangoni forces will be important in Skylab, but not on ground 
tests. 


C. CAPILLARY PUMPING IN TAPERED GAPS 

Upon melting the braze alloy in the M552 tapered gap specimens 
(SLN-4, MCN-4, MCN-5 and MCN-6) will have the configuration shown 
in Figure 8, Ignoring gravity effects and contact angle hysteresis 
r22 and 23], the pressure gradients determining flow direction are 
determined solely by simple capillary effects, as follows 1_24, 25 
and 26]: 

„ ^ 2a COS0 


2a COS0 



where the contact angle 0, the pressures P, P^ and P 2 and the radii 
of curvature R, and Ro are shown in Figure 8. The term a is 
the surface tension of the liquid vapor interface. 
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Sleeve 



Tube 


■ » 

Flow Direction 

IT 

(assuming 0 < 9 < ■^ ) 


FIGURE 8. 


TAPERED GAP CAPILLARY PUMPING 
FOR THE M552 EXPERIMENT 



By eliminating the common term P from these equations, the 
following results 

Pj - Pj = ZCT cose > 0 

because Rj is greater than R 2 . Thus the molten braze will always 
flow toward the narrow end of the tube. 

The same would be true when gravity is considered because 
ground tests were conducted vertically and g-levels aboard Skylab 
II were less than 10 Sg* i.e., the ratio of hydrostatic to capillary 
pressures are: 


< 3 X 10"^ 


Like gravity effects, contact angle hysteresis (advancing and receding 
angles differ) should also have negligible effects on final melt position, 
because hysteresis will only affect the rate of flow and not the flow 
direction. 

D. SURFACE TENSION MEASUREMENTS 

The values of surface tension, a, for both the Ag/Cu eutectic 
braze alloy BT720, and the same alloy with nominally 0.2% lithium 
addition, Lithobraze BT720, were measured as flSj: 



Temperature 

Surface Tension 

Braze BT720 

1173°K 

996 x 10“® N/cm* 


1273 

965 


1373 

956 


1473 

917 

Lithobraze BT720 

1273 

967 


1413 

942 


1473 

930 


Multiply by 10^ to get dynes/cm or ergs/cm^ 

The surface energies of the two Ag/Cu braze alloys are plotted 
as a function of temperature in Figure 9. With an error of 1% in 
a and 20°K in T, a least mean square line gives a slope of S =” do/dT 
= -0.25 X 10"= N/cm°K (-0.25 ergs/cm^OC). This can be compared 
with a value of S = -0.31 x 10"= N/cm°K for pure copper. 
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Curve ’653747-A 



Fig. 9 “Surface energy vs temperature Ag-Cu solder 
alloys. Al20^ substrate, argon atmosphere, [isj 
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E. FLIGHT RESULTS 

Radioisotope tracer mapping of Skylab samples SLN-2 and 
SLN-4 showed very interesting differences and similarities when 
compared with similar maps of ground-test samples. One significant 
difference was the enhanced braze alloy flow (with subsequent in- 
creased movement of tracer alloy) observed in the samples brazed in 
space under near-zero gravity conditions. Radioactivity was observed 
for both samples SLN-2 and SLN-4 well outside of the ferrule region 
on the double -scribed or tracer pellet end. With gravitational forces 
absent, capillarity of surface forces were \xnopposed, resulting in 
rather dramatic increased flow of the braze alloy. Another difference 
due to the absence of gravity was manifested in sample SLN-4, a taper 
annulus sample, which demonstrated the absence of a braze "leveling" 
effect noted on ground -test samples; islands of braze alloy filled re- 
gions having very large gap widths. In some respects, the radiation 
intensity maps for the Skylab samples SLN-2 and SLN-4 were very 
similar to those obtained from earth samples. The tracer alloy moved 
initially toward the thermally hot regions near the igniters. The 
tracer maps for all of NASA M-552 braze assemblies suggest that such 
tracing can provide a picture of the thermal history for any particular 
assembly as well as accurate braze alloy flow information. Since no 
other analytical technique would provide this vmique composite view 
of surface, gravitational, and thermal forces, radioisotope tracing 
proved useful in interpreting other metallurgical aspects of this ex- 
periment [2 7j. 

An observation of major importance is that true capillary flow 
can be developed in wide-gap apertures when gravity is absent. This 
had been predicted, and was a reason for inclusion of wide-gap speci- 
mens (SLS-3, 0.020 inch gap and SLN-4, tapering to 0.030 inch gap). 

There are practical as well as scientific implications to wide- 
gap capillary behavior. One principal consequence is that brazing, as 
a means for joining materials, is a much more versatile process in 
zero gravity space environment than it is on earth. Commercial ter- 
restrial capillary brazing is limited by gravity interference to gap 
clearances rarely exceeding 0.005 inch, which imposes serious con- 
straints on machining and fixturing. The absence of the need for 
close clearances and relatively smooth parallel surfaces, as required 
on earth, greatly extends the applicability of brazing as contrasted, 
for example, with welding, fastening, or other joining techniques, 
when structural fabrication is to be done in space environment. 

Some general characteristics of capillary flow were manifested 
in a review of metal distribution in all M552 Skylab specimens. The 
pressure in liquid near a small gap meniscus is smaller than near a 
wide gap meniscus; therefore pressure differences develop in complex 
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capillary systems tending to drive the liquid from wide gap to narrow 
gap locations. For example, when the braze gap was fairly narrow 
(SLS-1, 0.005 inch), the ring-grooves were almost drained of liquid; 
in a wide gap specimen (SLS-3, 0.020 inch), not only was the joint 
"starved" as mentioned above, but the wide gap was far less effec- 
tive in draining liquid from the ring-grooves (thus contributing in- 
directly to the "starvation" of the joint). In consequence, an 
observable, measurable meniscus developed in ring-grooves near 
wide gaps. For example, the upper ring-groove in specimen SLN-4 
is near a very wide (0.030 inch) gap, which was incapable of with- 
drawing much liquid from the groove and the residual liquid formed 
a meniscus according opportunity for observation of a radially sym- 
metrical liquid -vapor surface contour in the absence of gravity. This 
contour has been found numerically consistent with the following equa- 
tion which governs in the absence of gravity, provided the surface ten- 
sion is uniform from point to point on the meniscus surface [2s]. 



where 

Y = distance from the centerline of the specimen 
to a point on the meniscus, measured perpen- 
dicular to the centerline 

X = distance measured parallel to the centerline 
from the center of the ring-groove 

C = integration constant 

cr = surface tension 

AP = pressure difference prevailing across the 

meniscus surface (and constant from point to 
point along the meniscus surface) 

The conformance of the meniscus surface with the equation sup- 
ports the conclusion that the liquid -vapor surface tension is substan- 
tially uniform. This is of importance because the surface tension is 
quite sensitive to variations in temperature and surface composition. 

F. CONCLUSIONS 

The absence of gravity greatly extends the scope of brazing, and, 
thereby, the applicability of brazing to fabrication in space. In zero 
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gravity environment, the surface tension forces begetting capillary 
flow are unfettered, while on earth these forces must compete with 
gravity. Study of braze alloy distribution in Skylab specimens clearly 
indicates that dimensional tolerances, especially braze gap clearances, 
will be far less critical to joining operations in space than on earth. 

The practical significance of this fact, which had been predicted but 
never tested, can hardly be overemphasized. In space fabrication, 
many joints, which on earth would be produced by welding, should 
probably be brazed. 

Liquid-vapor boundary surfaces (menisci), and the flow of liquid 
metal under the impulsion of surface tension forces, are in close con- 
formance with what had been predicted for zero gravity environment. 
There were no unexplained effects, and, in the Skylab specimens, the 
surface tension of liquid silver- copper alloy appears to have been 
quite uniform. 
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SECTION 17. M553 SPHERE FORMING EXPERIMENT 


This chapter contains analyses of the physical forces for the M553 
experiment. The first section contains analyses of the driving forces, 
magnitude and pattern of natural convective fluid flows which occur in 
the Sphere Forming Experiment. Following sections address tra- 
jectories of released samples, shape dynamics, surface tension 
measurements, shrinkage pressure in freezing spheres and flight 
results. ° 

A. CONVECTION ANALYSIS 

As in the other M512 experiments, the physical forces and 
attendant fluid dynamics of the M553 experiment are important factors 
in determining the quality of the final product. Flow patterns in the 
molten material are important because all of the sample materials have 
low entropies of fusion [IJ. Thus, their solidification (microstruc- 
ture) is controlled by the rate of heat transfer removal [2] , which 
changes with the fluid flow [3 and 4] . The degree of flow will also 
determine the amount of mixing attained. If no or little flow were 
present, all heavier components would segregate to the lower portion 
of the spheres on earth, but to a negligible extent in the microgravity 
environment of space. Fluid flow can also affect the shape and re- 
lease of the specimen while it is retained on the ceramic holder. 

Application of dimensional analysis [5 and 6] to the governing 
equations for eb melting, coupled with ground-based and KC-135 ex- 
periments, should enable prediction of the extent of reduction or in- 
crease of motion in the molten metal and/or the change in flow pattern 
in electron beam melting in space. Possible physical forces which could 
induce fluid flow in the M553 experiment, and their causes, include: 

• Effect ive Gravity Force ; Resultant force on weld 
specirnen due to earth's gravity and centrifugal and 
Coriolis forces of orbiting spacecraft. 

• Lorentz Force ; Electromagnetic forces induced by 
passage of the electron beam current through the 
specimen. 

• Elect rostriction; Stresses induced when electrical 
permittivity changes with density. 
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• Magnetostriction : Stresses induced when permeability- 
changes with density 

• Electrostatic Force : Caused by presence of excess 
electrical charge (due to beam current and/or therm- 
ionic emission). 

• Surface Tension: Tangential stresses at vapor-liquid 
interfaces can be induced if surface tension depends 
on temperature and/or concentration. Surface tension 
will also cause pressure gradients across curved 
interfaces. 

• Density Differences Accompanying Phase Changes 

• Beam Force : Impinging electrons giving up their 
momentum. 

• Thermal Expansion ; Dilation and compression of 
fluids whose density changes appreciably with 
temperature can induce fluid flow. 

• V ib ration : Uncontrolled movement due to engine 
operation, astronaut motion, particle impacts, etc. 

• Centrifugal and Coriolis : Generated by disk rotation. 

• Vapor Pressure; Evaporating molecules impart 
momentum which leads to normal stresses at vapor- 
liquid interface. 

• Inertia Forces ; Tend to sustain induced motions. 

• Viscous Forces ; Tend to resist driving forces. 


Governing Equations 

The preceding forces, which could influence fluid flow and 
solidification, appear explicitly in the following conservation equa- 
tions which apply to formation of a molten pool by electron beam 
heating. 

Continuity (Mass Balance) : 

|^+ V ■ (pV) = 0 
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where 


p = density 

V = velocity vector 
t = time 

V = "grad" or "del" operator 


Momentum: 


p + p(v . V) V = -VP + ^ V^ V + p f 

o t 

Inertia Viscous Gravity 

Force Force Force 

+ J X ^ + Pe ^ 

Lorentz and Electro- Electrostatic 

magnetostriction Force 

Forces 

where 

P = pressure 
p = viscosity 
g = gravity vector 
J = electrical current density 
B = magnetic flux density 


p^ = excess charge density 
^ 

E = electric field density 
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Energy ; 


AT — ^ ? 

pCy + pCy (V . V) T = k T + $ + pQ 


where 


+ - {p t/k)v . y 

Thermal Expansion 
Source 


T = temperature 
Cy = heat capacity 
k = thermal conductivity 

= viscous dissipation function 
Q = internal heat sources 

= "conduction" current = o(E + V x B) 
o' = electrical conductivity 


P 

K 


_ 1 (l£- 

p 'at 


) 

P 


i ) 

p 'ap \ 


Maxwell's Equations : 


VXE = -3B/3t 
VXH = J + 3D/at 

V • B = 0 

V • D == 

where 

— i 

D = displacement current 
H = magnetic flux 



Constitutive Relations; 



B = ' H 

J = a' (E + V X B) + V 

p = p^ [^1 - pAT + KAp] 

whe re 

e ' = permittivity (dielectric constant) 

p ' = permeability 



Beam and Vapor Surface Tension Forces 

Pressure Forces 



where P^ ’ M ' . ^ are, respectively, the pres- 

sures, viscosities, and velocity components in Phases 1 and Z, a is 
the surface tension, R]^ and R2 principal radii of curvature of 

the surface, nj(i=l, Z, 3) are the components of the unit vector normal 
to the surface and directed into the interior of Phase 1, and s\umma- 
tion over a repeated index (k=l, Z, 3) is assumed, (See [l4j for details. ) 

Energy Balance, 

eo (T^ - ) - kn • V T = Q™ 

s 00 s EB 
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where 


e = emissivity 
<3 = Boltzmann's constant 

= surface temperature 
= environment temperature 

n = outward unit normal vector to free surface 
Q^B “ heat flux from electron beam 

Continuity of Velocity, 

n X (V*** - ) = 0 

f(l) = t. = 0 

Solid-Liquid Boundary Conditions ; 

Along stationary interfaces the no-slip condition (all components 
of velocity vanish) holds; whereas at moving interfaces (along melting 
or freezing fronts), a one -dimensional material balance yields 

Ui = (1 - Ps/Pl) d6/dt 

where 


m = fluid velocity normal to the interface 

6 = position of the interface relative to the 

origin of the spatial coordinate system 

= solid density 
Pb - liquid density 


This is the source of the density-difference-accompnahing-phase- 
change force for fluid flow. Also for a melting interface, an energy 
balance yields: 


-kV (T^ - Tg) = A f 
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whe re 


T = liquid temperature at the interface 
To = solid temperature at the interface 

O 

A = heat of fusion 


and 


The preceding equations are based on the following assumptions: 

• No influence of external fields on physical properties 

• No coupling between constitutive flux relations (e.g,, 
no Soret effect) 

• Single -component, Newtonian fluid 

• Constant physical properties 

• Bulk coefficient of viscosity vanishes. 

They are written in general vector form wherever possible, since they 
also apply directly to M551. It should be noted that only the differences 
in physical properties and geometry will differ in a dimensional analy- 
sis of M551 versus M553. 

Dimensional Analysis 

The controlling physical forces can be determined by nondimen- 
sioning each of the preceding equations and performing an order-of- 
magnitude conwarison on the various dimensionless groups which 
result [5 and 6J . The key to successful analysis is in choosing 
the proper reference values; i.e., since no freestream velocity exists, 
which forces do we equate to estimate a "characteristic" or "typical" 
velocity. Choosing the proper characteristic velocity is very im- 
portant, since the reference time, temperature, etc. usually depend 
on this velocity. 

Previous examinations of film showing both KC-135 flight ex- 
periments and groixnd tests for the M553 experiment indicated fluid 
velocities exceeded 300 cm/sec {p. B-5 of Reference 29). Continued 
M553 film analysis has since shown that flow velocities for nickel 
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specimens on both KC-135 and ground tests have been approximately 
^ possible forces affecting electron beam melting 

L those shown in Table 7 yield characteristic velocities 

of this order of magnitude. These driving forces consist of couplings 
tension, gravity or Lorentz (electromagnetic) forces 
with inertia forces. Previous studies based on faster velocities had 
indicated couplmg with viscous, rather than inertia, forces as con- 
trolling L^9J. This indicates the importance of choosing the correct 
characteristic velocity. Preliminary analysis also suggests that 
magnetostriction forces may also be important, but gross uncertainty 
in electromagnetic property data for liquid metals precludes any de- 
cision at the present time. 


Table 7 

PROBABLE CHARACTERISTIC VELOCITIES 
IN ELECTRON BEAM MELTING 



Velocity 

Controlling 
F orces 

Functional Form 

Value (cm/ sec) 
for Nickel in M553 

Inertia = Surface 
Tens ion 

(ST/pL)^/^ 

20 

Inertia = Lorentz 

(a' E^ L/p)^/^ 

1 

Inertia = Gravity 


18 at 1 
6 at 

Viscous = Lorentz 


7 


The dimensionless momentum equation which determines fluid 
flow in electron beam melting becomes (upon choosing the surface 
tens ion- inertia characteristic velocity). 
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N 


St 


9V 

9t 


N 


Oh 


(V ’ V ) V = - 


2N 


VP + V 


Oh 


(pg) + (VxB)xB 


N 


Oh 


Oh 


(IV. 1) 


where 


Ng^ = Stokes number = 


N 


= duration of process/residence time 
= Ohne s or ge number = ^ ^ ^ ^ "^o 

= viscous force/surface tension force 


N 


Bo 


z 

= Bond number = p^g^L /ST 


N 


MI 


=: gravity force/surface tension 
- Magnetic interaction number = 
= magnetic force/surface tension 


B^L 

o 


Zfx' ST 
^ o o 


Using physical property data for nickel and beam parameters 
for the M553 experiment, the equation reduces to 


10^ p ||- + 10^ p (V - V) V = 


- 1 0^ VP + V^ V 


+ G(pg ) + 10"^ (V X B ) X B 


where 

3 

G = 10 for earth gravity 

G = 10'^ for expected Skylab gravity. 

This order-of-magnitude analysis indicates that surface tension driven 
convection will occur both in ground tests and for Skylab conditions be~ 
cause Nqii « 1. This is confirmed by KC-135 M553 tests. Further- 
more, gravity driven convection will exist on ground tests, but will be 
negligible in the reduced gravity of Skylab. Thus different forces will 
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control convection on earth versus Skylab. The preceding analysis 
also indicates that electromagnetic or Lorentz forces will be negligible 
with regard to causing fluid motion. 

The surfaV:e tension driving force considered above is actually a 
surface tension gradient caused by radial and lateral temperature 
gradients. Ignoring convection, gradients of at least several hundred 
deg^rees Celsius per centimeter have been predicted during melting 
[ 31 J. It can be shown by dimensional analysis (equate inertial and 
viscous terms) that the motion caused by the initial temperature 
gradients will occur in less than 0.1 second and can persist for 60 
seconds after removing the driving force. This means that there will 
be some fluid motion during solidification if the M553 specimens freeze 
after 30 to 40 seconds as predicted. 

Values of the pertinent dimensionless groups for the remaining 
M553 materials are given in Table 8- -As can be seen, no significant 
changes are evident from the preceding conclusions for pure nickel. 


Table 8 

M553 DIMENSIONAL ANALYSIS 


Material 

^Oh 

Bo 

1 

* 

^Bo 

N t 

Bo 

'^Oh 

^Oh 

^Oh 

Ni-Sn 

4.4 X 10'^ 

1.50 

2292 

3437 

.34 

Ni-Ag 

5.0 X 10"'^ 

1.53 

1991 

3043 

.30 


.4 





Ni-Cu 

4.45 X 1 0 

1.54 

2107 

3241 

.32 


With Earth gravity (9.8 m/sec^) 
^With 10 Earth gravity 


It should also be mentioned that the effects on fluid motion of sur- 
face free charges (excess electrostatic charge) and related electric and 
magnetic forces at the drop surfaces have been examined. The results 
of the examination indicate that these surface charges are negligible 
d riving forces for fluid motion. These results were also reviewed by 
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Professor J. R. Melcher of MIT's Electrical Engineering Department, 
and he was in full agreement with these results"'. Furthermore, 
from the manner in which patches of surface contaminants moved 
about in the ground films, it is apparent that surface tension driven 
flows due to concentration gradients are also important in this ex- 
periment. No reliable data exist, however, on what these impurities 
are, nor what the value of surface tension gradient with composition 
is. Thus the Marangoni effects are limited to thermal differences 
only in this study. 

B. M553 TRAJECTORY ANALYSIS 

One of the most important aspects of the M553 experiment is the 
motion of the sphere after deployment since this determines what free- 
float time will exist. In order to predict the motion, a trajectory pro- 
gram has been developed that includes allowances for: 

• Electron beam force 

• Deployment velocity (spring) 

• Vaporization force based on 3-D temperature 
history 

• Sky lab orbit considerations 

• M512 position considerations, and 

• Allowances for additional forces as independent 
subroutines (KC-135 g-level data). 

In general a solution to the equation of motion 



i (t) 


is obtained numerically via Simpson's rule subject to the following 
boundary conditions: 


R 


0 and 


dR 

dt 


V at t = 0 
o 


Integration of the equation of motion twice yields 

t r) 

R = t v^ + jy 7" (e ) de dn 

o o 


Private Communication, Huntsville, Ala., 
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1 1 June 1973. 



which produces the desired trajectory in chamber coordinates centered 
on the specimen. This computed trajectory is then transformed to 
camera coordinates to show the path as seen in the film, i.e., 


R'(t) = [a] s (t) 

I 

where R (t) is the projected trajectory and Fa] is the transformation 
matrix obtained from the hardware. For cani^ra 1 


r 


A = 

.932 

-.203 

-.363 

-.521 

0 

-.829 


.301 

-773 

-.560 

and for camera 2 

A 

^ .661 
-.224 

-.183 

-.974 

-.728 

.041 


-.716 

,136 

.685 


force was discussed in Reference 31 and Table 9 gives 
the magnitude of this force for two eb currents. To establish the 
motion as a result of the eb force, the conversion efficiency and eb 
cutoff time must be known. This time corresponds to the time delay 
associated with automatic shutdown of the eb due to sting melting 
Original estimates were from 10 to 20 milliseconds; however, high 
speed movies (1000 frame/sec) taken during the M553 ground-based 
test at MSFC on 20 October 1972 show conclusively that longer times 
exist. During their test the spring loads were varied on the shutdown 
mechanism to experimentally establish the time interval the eb is on 
after melting occurs. It was found that for a 150-g spring load a cut- 
off time of 410 milliseconds was measured while 220 milliseconds were 
observed for a 200~g load. This cutoff time translates into a free- 
float time of approximately 12 seconds for a 150-g load and 23 seconds 
for the 200- g load. Thus the expected free-float time is drastically 
reduced from the 250 seconds corresponding to the original 20 milli- 
seconds estimate. ® 


An analysis was performed to determine if the cohesive energy- 
generated kinetic energy is sufficient to cause deployment in the 
absence of external forces. This was done to allow initial conditions 
to be calculated. To reduce the complexity of the problem, the fol- 
lowing assumptions were made; 
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Table 9- 

TOTAL FORCE DUE TO ELECTRON BEAM 


F 

e 


= I 



V 




m 


F 

F 

I 

e 

e 

V 

e 

e 

(amp) 

(coulomb) 

(kg) 

(V) 

(N) 

(dyne) 

0.05 

1.6(10)'^^ 


2(10)^ 

2.38(10)"^ 

2.38 

o.io 

1.6(10)“^*^ 

0.9(10)"^° 

2(10)^^ 

4.76 (10)“^ 

4.76 


• Zero gravity field 

« No sting melting 

• Sting pulled down before breaking 

• Perfectly spherical free surface, and 

• Contact angle of 120 degrees. 

An energy balance was made between the initial and final states 
and from this the kinetic energy is found as 

KE = (a^/> + 

a I A - A 

LS ^'^LS LS ) 

This kinetic energy was calculated as a function of penetration 
(H) of liquid into orifice. The value of H simulates the effect of the 
spring in pulling the liquid down via adhesion between liquid and solid 
sting. The results of this calculation are 
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Penetration, H 
(cm) 

0.0 

0.1 

0.5 

1.0 


Kinetic Energy 
(ergs) 

0.0 

62.7 

293.0 

597.0 


The value of the kinetic energy must then be greater than the 
work of adhesion (or energy of adhesion) for deployment to occur. The 
work of adhesion was found to be 265 ergs; thus for depressions of more 
than 0.5 cm deployment will occur due to cohesive effects alone. These 
criteria are included in the trajectory computer program to allow initial 
conditions for the numerical integration to be determined. 

Tables 10 and 11 show the surface temperature distribution at 

3.0 and 3.5 seconds, respectively, for Ni and Ni-12% Sn. These data 
were used in conjunction with the Langmuir theory discussed in 
Reference 31 to calculate the net vaporization force as a fcuiction of 
time. A sample calculation is presented in Table 12 for both Ni and 
Ni-12% SN for the first second after eb cutoff. It is seen that between 
3.5 and 4.25 seconds the vaporization force is of the same order of 
magnitude as the eb momentum force. During the first half second 
the calculations yield a large vaporization force (comparable to surface 
tension forces) and represent an intolerable situation as far as the free- 
float time is concerned. However, due to the large temperature gra- 
dients existing, the resulting surface tension driven convection will 
tend to reduce these values. 

The Skylab orbit effects are included in the analysis by the method 
of Parker and Gatewood [32j. The computer program presented in 
Reference 31 was used as a subroutine in the present analysis to couple 
the effects of M512 position and orbit to the physical phenomena occur- 
ring in the vacuum chamber. 

C. CONCLUSIONS BASED ON TRAJECTORY COMPUTATIONS 

It has been shown that free-float times from 32 to 48 seconds 
(depending on the material) are required for complete containerless 
solidification to occur. The magnitude of each of the physical forces 
was discussed and general comments made concerning the effect of 
each on free-float time in the previous section. This subsection sum- 
marizes some of these conclusions. 
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TABLE 10 


TEMPERATURE DISTRIBUTION AT 3.0 SECONDS 



Nickel 

Nickel - 12% Tin 

Node 

No, 

Temp. 

(°C) 

Vapor Pressure 
(mm Hg/ cm^) 

Temp. 

(°C) 

Vapor Pressure 
(mm Hg/ cm^) 

36 

37 

2036 

1463 

.3959 

.0006 

2013 

1537 

.7001 

.0069 

38 

2193 

.464 

2157 

.6965 

39 

1900 

.0395 

1903 

.0962 

40 

1728 

.0067 

1757 

.02 53 

41 

i 

1677 

.0037 

1716 

.0168 

4Z 

1677 

.0037 

1716 

.0168 

43 

1728 

.0067 

1757 

.02 54 

. 44 

1900 

.0400 

1903 ! 

.0963 

45 

2193 

.4644 

2157 

.6965 

46 

1801 

.0202 

1815 

.0602 

47 

1695 

.0063 

1731 

.0266 

48 

1609 

.0022 

1663 

.0132 

49 

1573 

.0014 

1634 

.0096 

50 

1573 

! .0014 

1634 

.0096 

51 

1609 

1 .0022 

1663 

.0132 

52 

1695 

.0063 

1731 

.0266 

53 

1801 

' .0201 

1815 

.0602 

54 

1633 

.0030 

1680 

.0159 

55 

1584 

.0016 

1642 

.0105 

56 

1537 

.0009 

1605 

.0069 

57 

1513 

.0006 

1586 

.0056 

58 

1513 

.0006 

1586 

.0056 

59 

j 1537 

.0009 

1605 

.0069 

60 

1584 

.0016 

1642 

.0105 

61 

1633 

.0020 

1680 

.0158 

62 

1535 

.0006 

1602 

.0049 

63 

1515 

.0005 

1585 

,0041 

64 

1492 

.0003 

1567 

.0033 

65 

1480 

.0003 

1557 

.0030 

66 

1480 

.0003 

1557 

,0030 

67 

1492 

.0003 

1567 

.0033 

68 

1515 

.0005 

1585 

.0041 

69 

1535 

.0006 

1602 

.0049 











TABLE 11 

temperature distribution at 3.5 SECONDS 


Node 

No. 

Nickel 

Nickel -12% Tin 

T emp. 
(°C) 

Vapor Pressure 
(mm Hg/cm^) 

Temp. 

(°C) 

Vapor Pressure 
(mm Hg/cm2) 

36 

1689 

.01249 

1678 

.0332 

37 

1601 

.00429 

1624 

.0185 

38 

1685 

.0041 

1676 

.0111 

39 

1679 

.0038 

1674 

.0109 

40 

1672 

.0035 

1671 

.0105 

41 

1667 

.0033 

1670 

.0103 

4Z 

1667 

.0033 

1670 

.0103 

43 

1672 

.0035 

1671 

.0105 

44 

1679 

.0038 

1674 

.0109 

45 

1685 

.0041 

1676 

.0111 

46 

1675 

.0050 

1672 

.0144 

47 

1668 

.0046 

1 669 

.0140 

48 

1658 

.0041 

1665 

.0135 

49 

1652 

.0038 

1662 

.0131 

50 

1652 

.0038 

1662 

.0131 

51 

1658 

,0041 

1 665 

.0135 

52 

1668 

.0046 

1669 

.0140 

53 

1675 

.0050 

1672 

.0144 

54 

1657 

.0040 

1663 

.0131 

55 

1651 

.0037 

1660 

.0128 

56 

1642 

.0033 

1657 

.0123 

57 

1636 

.0031 

1654 

.0200 

58 

1636 

.0031 

1654 

.0200 

59 

1642 

.0033 

1657 

.0123 

60 

1651 

.0037 

1660 

.0128 

61 

1657 

.0040 

1663 

.0131 

62 

1634 

.0022 

1649 

.0083 

63 

1630 

,0021 

1648 

.0082 

64 

1624 

,0020 

1645 

.0080 

65 

1621 

.0019 

1644 

.0078 

66 

1621 

.0019 

1644 

.0078 

67 

1624 

.0020 

1645 

.0080 

68 

1630 

.0021 

1648 

.0082 

69 

1634 

.0022 



1649 

.0083 
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Table 12 

SUMMARY OF VAPORIZATION FORCE CALCULATION 


Nickel 

Vaporization 

Force 

Time (sec) 

(dyne s) 

3.00 

3.25 

3.50 

3.75 

4.00 

4.25 

£ 

X 

882 

55 

8 

1.5 

0.3 

.07 

f 

y 

1480 

126 

25 

5.8 

1.3 

.3 

f 

z 

292 

25 

5 

1.0 

0.25 

.06 

Nickel- 12% Tin 

f 

X 

1400 

102 

9 

1 

0.16 

.02 

f 

y 

2470 

240 

41 

6 

1.0 

.16 

f 

z 

485 

47 

8 

1 

m 

.03 


• Electron beam force for a 460-g spring load has a dura- 
tion of approximately 100 milliseconds and gives rise to 
a float time of 40 seconds. 

• Kinetic energy of sphere due to spring effects (deploy- 
ment technique) gives rise to float times of greater 
than 50 seconds. 
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• Vaporization force produces float times that are probably 
the critical link in process. Best estimates of the cor- 
responding float time between l/4 second after eb cutoff 
and occurrence of an isothermal outer layer range from 

5 to 20 seconds. For some cases the time could even be 
considerably less than 5 seconds. 

• Skylab orbit effects give rise to 29 to 46 seconds of 
float time. 

• Superposition of ail of these forces produces an additive 
(as opposed to canceling) effect in most cases. Trajectory 
calculations show that even with conservative vaporization 
forces float times of less than 15 seconds will occur. 


Lockheed documented in December 1972 [33j an early projec- 
tion of anticipated float times less than the required containerless 
solidification times and again in April via the Science & Engineering 
Information Technical Note [34]. In the latter report a conservative 
estimate of less than 20 seconds w^s forecast. 

Another potential problem area, discussed in the 27 March 1973 
presentation at MSFC, is that even with the relatively large forces on 
the sphere during the critical time of deployment, the possibility of 
"sticking" exists. Films of the KC-135 flight test showed that an 
adherence to the ceramic occurs on occasion that finally results in 
deployment due to acceleration transients in the ballistics trajectory 
of the aircraft. A typical sequence is shown in Figure 10. The events 
are referenced to the time corresponding to what appeared to be the 
fully molten condition, t=0. For the ensuing 0.300-second, that adher- 
ence to the ceramic occurred with periodic motion due to the accelera- 
tion transients is evident in the film sequence. At 0,310-second, 
deployment occurred with the trajectory shown in Figure 10. The 
important point remains, however, that without the acceleration 
transients (relatively large compared to the Skylab gravitational en- 
vironment) specimen retention could have occurred. This adhesion is 
probably due to contamination of the ceramic from previous melts. 

The conclusion is that the possibility of unintentional retained speci- 
men may occur in the actual Skylab mission. 

D. SHAPE DYNAMICS 

This subsection presents a summary of calculations made for 
Ni and Ni-12% SN based on the thoughts outlined in Reference 31. 

First, the formation times are shown in Table 13 where it is seen 
that the transformation to minimum energy occurs very rapidly. How- 
ever, due to inertial considerations, the spherical free surface is 
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t = 0 (sec) 


t = .250 (sec) 




t = .275 (sec) 


t = .290 (sec) 




t = .310 (sec) 


FIGURE 10. TRAJECTORY SEQUENCE - MARCH 1972 KC-135 FLIGHT 
(NICKEL) (CONTINUED) 
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t = .325 (sec) 


t = .360 (sec) 



FIGURE 10. TRAJECTORY SEQUENCE (CONCLUDED) 
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"overshot" and oscillations are set up which continue until dissipated 
by viscous damping. Table 15 gives the oscillation frequencies while 
Table 14 shows the time required for viscous effects to dampen an 
initial perturbation to 0.01% and 1% of its initial value. 


Table 13 

FORMATION TIMES FOR M553 MATERIALS 


M553 Materials 

P 


R 

T 


(g/cm^) 

{dyne/ cm) 

(cm) 

(sec) 

Ni 

7.85 

2050 

.3175 

.011 

Ni-12% Sn 

7.80 

1500 

.3175 

.013 


Table 14 

FREQUENCY FOR M553 MATERIALS 


Material 

Frequency 


(Hz) 

Ni 

41 

Ni-12% SN 

35 


Table 15 


DECAY TIME FOR M553 MATERIALS 


Material 

P 

P 

R 

Time 

Time 

Time 


(g/ cm-sec) 

(g/cm^) 

(cm) 

(constant) 

(.01%) 

(1%) 

Ni 

.05 

7.85 

.3175 

3.14 

28.8 

14.5 

Ni-12% Sn 


.043 

7.80 

.3175 

3.62 

33.0 

16.7 
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E. SURFACE TENSION MEASUREMENTS 

The values of cCj^yi and W. determined from sessile drop 
experiments ^1 5j are given in TableU6. 


Table 16 

WORK OF ADHESION (W^) FOR NICKEL AND ITS ALLOYS 

ON Af 203 (x 10-5 N/cm) 



Temp. (°K) 

*^LV 

9 (deg) 

^A 

Nickel 

1773 • 

1719 

135 

503 


1973 

1676 

137 

450 


2068 

1659 

(113) 

(1011) 

Ni-1% Silver 

1793 

163 7 

135 

479 


1983 

1586 

136 

445 


2138 

1564 

128 

601 

Ni-12% Tin 

1683 

855 

130 

305 


1983 

806 

133 

256 


2093 

782 

130 

279 

Ni-30% Copper 

1743 

1506 

126 

621 


2003 

1468 

130 

524 


2073 

1411 

133 

449 


F. SHRINKAGE FORCES 

When spherical melts solidify, the innermost melt is the last 
to freeze. The advancing outer solid shell leads to hydrostatic ten- 
sions, even negative pressures, in the residual liquid and the general 
plastic collapse of the spherical casting. The contraction (about 3% 
for most metals) is caused by solidification shrinkage. Several recent 
analyses have been performed for this configuration [35 through 
37], and two extreme cases which could exist are shown in Figures 
11 and 12, In Figure 11 the residual liquid is considered to have 
undergone cavitation. The homogeneous fracture pressure of most 
liquid metals is thought to be on the order of -10'^ atm; however, 
heterogeneous nucleation may be possible at less negative pressures. 

The tensile pressures within the residual liquid can be calcu- 
lated assximing various growth rates and outer shell deformation modes: 
elastic-plastic, Bingham, viscous or creep flow. Using the creep 
flow model, negative pressures in the liquid core of 1 cm radius nickel 
sphere are shown in Figure 13 using initial pressure of 1 atm. Nickel 
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FIGURE 11 . UNSTRAINED OUTER CASTING AND 
HOMOGENEOUS NUCLEATION OF 
RESIDUAL LIQUID UNDER EXTREME 
NEGATIVE PRESSURES [35] 



FIGURE 12 . STATE OF A SOLIDIFYING SPHERICAL 
CASTING with EXTERIOR SHELL 
DEFORMATION [SSj 
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FIGURE 13. 


HYDROSTATIC TENSIONS IN THE LIQUID CORE 
OF 1 cm RADIUS CASTINGS OF ALUMINUM, 
COPPER, NICKEL. AND IRON CALCULATED 
FROM THE CREEP MODEL [36J 


spheres in the M553 experiment are 0.32 cm in radius. The tensile 
pressures will be lesser in wide freezing range alloys 

Thus the value of possible shrinkage pressure developed withm 
spherical solidifying castings has been evaluated for various materials 
exhibiting various deformation characteristics. The maximum at- 
tainable stresses in aluminum, copper, nickel and iron are in the 
range -10'^ to -10“^ atm. The exact values are uncertain because 
the solidification rate of a sphere is not accurately known, but it is 
probable that homogeneous nucleation is not a feasible mechanism 
for the initiation of porosity in solidifying metals. These extremely 
low pressures may help explain the unexpected results being found by 
other M553 contractor investigators. 


57 



G. FLIGHT RESULTS 


A compilation of fluid velocities (on the molten sphere surfaces) 
measured from movies taken during Skylab, KC-135 aircraft and ter- 
restrial processing is shown in Table 17, These velocities are meas- 
ured by tracking oxide or shiny contaminant patches as they moved over 
the molten surfaces both before and after complete melting. Several 
KC-135 measurements were obtained even after release from the 
ceramic pedastals. Shown in the table are the surface velocity, 
processing environment, gravity level, date of processing, filming 
information, specimen identification and time of measurement rela- 
tive to the melting sequence. Beam power data for all measurements 
were approximately 50 mA and 20 kV. The velocity measurements are 
accurate to within 20%. The variations in results are probably due to 
differences in material compositions and beam settings (i.e., focus, 
impact point, voltage and amperage) rather than gravity effects. 

As shown earlier only surface tension and gravity forces are 
predicted to produce these velocities. The one velocity measure- 
ment obtainable from the Skylab movie was an order of magnitude 
lower than that obtained from some low gravity KC-135 aircraft and 
ground-based movies, but was similar to ground and KC-135 measure- 
ments also obtained long after complete melting. Thus it can be con- 
cluded that no significant difference in magnitude of convection occurred 
between Skylab, KC-135 and ground processing, which reinforces the 
theory of surface tension as the dominating force for fluid motion. 

In addition to direct movie evidence of vigorous melt stirring, 
several other effects of convection were evident in metallurgical 
evaluations of Skylab specimens conducted by other investigators. The 
grain multiplication near the large shrinkage pore in the cap region of 
Skylab sample SL-1.6 (nickel-tin material) indicates vigorous flow 
probably due to solidification shrinkage [38 ]. Alumina was spread 
over the surface of SL-2.?[ 39 ] a release -type nickel-tin specimen, 
which supports the existence of vigorous convective stirring to trans- 
port the alumina from the pedestal. Similar findings of alumina trans- 
port were found on SL-1.8[ 40 ], a nickel-copper specimen. 

Several Skylab samples possessed a large spherical or oblong 
cavity. These large interior cavities were usually associated with a 
"cap region" opposite the flat, epitaxial side of the specimens. These 
large pores are no doubt produced by the large shrinkage forces pre- 
dicted to exist during Skylab processing. 
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Table 17 

M553 SURFACE VELOCITY MEASUREMENTS 


Processing 

Environment 

Gravity 

Level 

(g) 

Date 

Film 

j. 

(fps) 

Sphere 

Material 

Fluid 

Velocity 

(m/sec) 

Time of 
Measurement 

Skylab 

o 

1 

6/15/73 

Color, 

Ni-Ag 

0.02 

Very long after 




24 



melted 

Ground 

1 

4/30/73 

B&W, 

Ni-Ag 

0.35 to 

From early to 




3000 


0.40 

complete melt 

KC-135 

to 3 

7/25/72 

Color, 

Ni-Ag 

0.20 

Before full melt 


10 


200 


0.10 

Right after full 







melt 






0.06 

Long after full 







melt 

Ground 

1 

10/20/72 

B&W, 

Ni 

0.15 

Before 




1000 


0.20 

complete 






0.33 

melt 

G round 

1 

1972 

Color, 

Ni 

0.12 

Half melted 




200 





y 2 

Ig = 9.8 m/sec 

fps = frames per second 
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H. CONCLUSIONS 


The primary conclusion of the preceding discussions is that 
significant flow and stirring in the molten spheres (during and after 
eb melting) occurred in both terrestrial and space processing as pre- 
dicted. The magnitude of this convective fluid motion was the same 
in both environments. The flow patterns, however, could not be 
quantitatively compared because the filming rate was too slow. Also, 
altered microstructure, increased grain refinement and the appearance 
of a single, large interior shrinkage pore have been found in the Skylab 
specimens as predicted. These effects can be traced to a longer dura- 
tion containerless freezing than in ground and KC-135 flights; i.e., the 
surfaces opposite the expitaxial nucleation site were cooler than in the 
non-Skylab specimens (assuming release). These surfaces nucleated 
much sooner than their non-Skylab counterparts, thereby entrapping a 
larger enclosed pool of residual liquid. The solidification shrinkage 
forces, which cause fluid flow and tensile stresses (increased shrinkage 
porosity), were therefore much greater in the Skylab specimens even 
though all of the samples appeared to have some degree of unplanned, 
heterogeneous nucleation (retainment of unmelted material or solidified 
prior to release). 

It should be noted that a repeat of this experiment in the micro- 
gravity environment of space as is, except for longer, higher speed 
camera coverage, would yield much more hydrodynamic information. 
This information is important because hydrodynamics controls the re- 
lease of the specimens and affects the solidification and microstructure 
properties. The release mechanism is of utmost importance because 
it will control the duration of free float time in the absence of acoustic 
or electromagnetic levitation. It is also recommended that a new ex- 
periment be performed with no active mechanical forces (springs) for 
release, i.e., allow only the vaporization and electromagnetic forces 
to initiate release. The slight forces could be counteracted by the 
small residual gravity field to yield very slow release acceleration and 
extremely long free -float times. 
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